Dynamical effects of an unconventional current-phase relation in YBCO dc-SQUIDs 
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The predominant d-wave pairing symmetry in high temperature superconductors allows for a va- 
riety of current-phase relations in Josephson junctions, which is to a certain degree fabrication con- 
trolled. In this letter we report on direct experimental observations of the effects of a non-sinusoidal 
current-phase dependence in YBCO dc-SQUIDs, which agree with the theoretical description of the 
system. 
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It is well established |l| that the wave function of a 
Cooper pair in most cuprate high-temperature supercon- 
ductors (HTS) has a d-wave symmetry. Its qualitative 
distinction from e.g. the anisotropic s-wave case is that 
the order parameter changes sign in certain directions, 
which can be interpreted as an intrinsic difference in the 
superconducting phase between the lobes equal to tt. 

The latter leads to a plethora of effects, like formation 
of Andreev bound states at surfaces and interfaces in cer- 
tain crystallographic orientations 0, S The current- 
phase dependence Is(<f>) m Josephson junctions formed 
by (id-junctions, as well as by srf-j unctions comprised of 
a cuprate and a conventional superconductor, depends 
both on the spatial orientation of the d-wave order pa- 
rameter with respect to the interface, and on the qual- 
ity of the latter 0, 0, S • Time- reversal symmetry 
can also be spontaneously violated and thus spontaneous 
currents generated 

[3 111 113 Another effect can be 
doubling of the Josephson frequency |6J, IjJjl ll4 • 

In this letter we report on experimental observations 
of strong effects of an unconventional current-phase rela- 
tion on the dynamics of two (id-junctions integrated into 
a superconducting interference device (SQUID) configu- 
ration. 

Since /<?(</>) must be a 27r-periodic odd function, it can 
be expanded in a Fourier series. In the most cases only 
the first two harmonics give a significant contribution to 
the current: 



l[ sin (f> - ll 1 sin 2<f> 



(1) 



In Josephson systems of conventional superconductors 
the second harmonic will usually be negligible [lfij but 
in (id-junctions the second harmonic may dominate. If 
PJ > the equilibrium state is no longer = but 
becomes double degenerate at (j> = ± a,TCCos(l[ /2I 11 ) — > 
7r/2. The system can then spontaneously break time re- 
versal symmetry by choosing either state. Spontaneous 
currents as well as fluxes can be generated in this state. 



The potential will have the shape of a double well and 
there are reasons to believe that it will be possible to 
observe quantum-coherence in this system. The presence 
of a 2 nd harmonic in the current-phase relation(CPR) of 
a dd-junction was confirmed by IPichev et al. [8j- 

A non-sinusoidal CPR of the junctions will change the 
dynamics of a dc- SQUID 16] . Regarding the junctions as 
magnetically small, the supercurrent through the SQUID 
in the presence of an external flux $> x = $o ■ (</ , s/27r) can 
be written as 



I s {4>, <t> x ) = l' cX sin0 - l[{ sin(20) + Z C J . 
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The critical current through the SQUID is given by the 
usual expression I c {<f> x ) — max.* I s {4>i 4>x)- The time- 
averaged voltage over the SQUID in the resistive regime 
is readily obtained in the resistively shunted junction 
(RSJ) approximation. By introducing 5jd, <p x ] = 4>2~ 4>i 
and applying the same method as in |17| with the nec- 
essary generalizations, we obtain for the average voltage 
over the SQUID: 



v- L = 



Gi + G 2 
2tt 



-h 



2e at 



(3) 



Here are the normal conductances of the junc- 
tions, and 

5 + ^ + ?( /2 ^ - 5 /2) - h(4> + 5/2)) = 0, (4) 

gives the difference, 6, in phase drops across each junc- 
tion. In deriving il-SI 111 we have assumed that the induc- 
tance L is equally divided between the SQUID arms. We 
have also neglected the spontaneous magnetic fluxes in 
the dd-junctions, due to their small amplitude [ill Ilij. 



2 




-2 -1.5 -1 -0.5 0.5 1 1.5 2 



T x 

FIG. 1: The results of simulations of the I c — (j> x and V — <j> x 
dependence for a dc-SQUID with f ol = 1, I* 2 = 0.1, I'J = 0.2 
and 1^2 = 0.4 (arb. units). The different curves correspond 
to bias currents in the range / = I ! cl to I = 5/^1 -We assume 
L = and Gi = G 2 . 

Though (0} is only explicitly solvable in the limit L — > 0, 
it always yields £[— 0, — <?y = — <5 [0, <?y . This means that 
the usual inversion symmetry is retained. 

The results of numerical calculations based on J3J and 
@ are shown in fig^ The cusps in the critical current 
correspond to the points at which the global maximum 
in switches from one local maximum to another |l6| . 
Note the quasi-$o/2-periodicity of the current isolincs 
in the V — cf> x picture, reflecting the current-phase de- 
pendence and their shift along the "I^-axis, which 
depends on the sign of the bias current (as it must to 
maintain the central symmetry with respect to the ori- 
gin) . The shift does not depend on the magnitude of the 
current since we neglect the self-inductance. For large 
biases the $o _ periodicity is restored. Indeed, as the bias 
grows, one set of minima of the washboard potential, 
U = (h/2e)[-I z coscf) + {I 11 1 2) cos 20 - I(f>], disappears 
first unless the first harmonic I 1 is exactly zero. 

We have fabricated and studied a large number of dc- 
SQUIDs. The samples were fabricated from 250 nm thick 
YBCO-films deposited on SrTiC>3-bicrystals. The grain- 
boundary junctions (GBJs) are of the asymmetric [001]- 
tilt type with the misorientation angle of 45° (0° — 45° 
GBJ). For more information on GBJs see for example 
reference [lf| . 

The pattern was defined using E-Beam lithography 
and then transferred to a carbon mask employing a multi- 
step process. Finally, the YBCO is etched through the 
mask using ion-milling. This scheme allows us to fabri- 
cate high-quality bicrystal junctions as narrow as 0.2 /im, 



as has been reported elsewhere |20j. In the SQUIDs un- 
der investigation the junctions are nominally 2 fim wide; 
hence the fabrication-induced damage of the junctions is 
small. 

The measurements were done in an EMC-protected en- 
vironment using a magnetically shielded LHe-cryostat. 
However, the magnetic shielding is imperfect, as is evi- 
dent from the fact, that the expected zero-field response 
of our SQUIDs is not exactly at zero. The measuring elec- 
tronics is carefully filtered and battery-powered whenever 
possible. In order to measure the dependence of the crit- 
ical current on the applied field we used a voltage dis- 
criminator combined with a sample- and- hold circuit. All 
measurements reported here were performed at 4.2K. 

The SQUID loops are 15x15 /zm 2 . The numerically 
calculated inductance [2lJ is approximately 25 pH, yield- 
ing the factor (3 = 2irLI c /$o between 0.5-2. 

The SQUIDs were largerly non-hysteretic with a re- 
sistance of about 2 f2. The measured critical current 
varies from sample to sample but is in the range of tens 
of microamperes giving a current density of the order of 
J c = 10 3 A/cm 2 . The estimated Josephson penetration 
length A ,7 = ^o/V^^ajc^L is approximately 2 /im in 
all junctions, which means that the junctions are mag- 
netically short. This is supported by the quasi-period 
of the pattern in fig0 being close to the expected value 
0o/2At,w 17]. The differential conductance curves do 
not show any trace of a zero bias anomaly (ZBA), as is 
expected for 0° — 45° GBJs. ZBAs has been observed by 
other groups in GBJs with other orientations 0. 

The critical current is plotted as a function of applied 
magnetic field for two SQUIDs in fig|3 The result is 
in qualitative agreement with theory if we assume that 
the SQUID junctions have different ratios of the 1 st and 
2 nd harmonics of the critical current. This assumption 
is supported by the fairly small modulation depth (it is 
easy to see from equation © that I c would go exactly to 
zero in a SQUID with junctions of identical I C 2/Ici)- 

We can fit the data to equation (0), if we compensate 
for the residual background magnetic field and assume 
that we have a small excess current (of the order of a 
few /iA) in the junctions. The fitting parameters again 
confirm that there is a large asymmetry between the arms 
of the SQUIDs. Note, that the model does not consider 
the flux penetration into the junctions, 

The result for fields of the order of mT is presented in 
fig'Gl which shows the 7 c -modulation of the SQUID en- 
veloped by an anomalous Fraunhofer-pattern quite sim- 
ilar to what has been reported by other groups [22I E^] 
for 0° — 45° GBJs. Note the inversion symmetry of the 
pattern with respect to the origin. That the global max- 
imum is not in the center can be explained in several 
ways; it has been shown for example that this could be 
due to the presence of so-called 7r-loops in the junction 
interface [24| . 

Figure 01 shows the V — independence of one of the 
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FIG. 2: Critical current as a function of magnetic field at 
4.2K. The dashed box indicates the area plotted in fig. [3Ji. 
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FIG. 3: Critical current as a function of applied magnetic 
field for two different SQUIDs that are nominally identical. 
The solid line represents the fitted expression. The fitting 
parameters are as follows: (a) i^i = 9 [J.A, 1*% = 0.3 fiA, I*{ — 
3.7 fiA and i" = 22.7 fiA (b)/ii = 7.8 fiA, l' c2 = 3.0 fiA, 
ill = 5.3 fiA and ill = 4.3 /iA . In both cases the fit has 
been adjusted with respect to the residual background field 
and the excess current of the junctions. 



FIG. 4: Voltage modulation as a function of applied magnetic 
field for the SQUID whose I c — B is shown in fig. [3^. The 
pattern is again inversion symmetric. Note the sign change 
at 100 fiV which we believe is due to a LC-resonance in the 
SQUID loop. 



SQUIDs. The pattern is again field inversion-symmetric. 
The overall structure is the same as in the model depen- 
dence of figd but there is also an additional shift due 
to self-field effects, which depends on the magnitude of 
the bias current and corresponds (at maximum) to a flux 
~ 0.1$o- In a beautiful experiment a similar dependence 
was recently observed by Baselmans et al in a Nb-Ag- 
Nb SNS junction where current-injectors were used to 
change the occupation of current-carrying states in the 
normal region [25| . A deviation from the model occurs at 
y=100 //V where the minima and maxima switch. This 
is probably due to an LC-resonance in the SQUID. Tak- 
ing L — 25 pH, this would require C — 0.8 pF, which 
agrees with our measurements on single junctions 

Remarkably, the observed offset of the V — B- 
characteristics with respect to the direction of the bias 
current appears to be a much more robust manifestation 
of the presence of a second harmonic of the Josephson 
current, than the shape of the I c — B curves itself. We 
observed the shift even in SQUIDs with the smallest junc- 
tions down to 0.5 fxm wide, where the deviations from the 
usual sinusoidal CPR were not obvious from the I c — B 
dependence. 

Generally the nature of the transport through a GBJ 
will depend on its transmissivity D. Il'ichev et al. @ 
have reported values of D as high as 0.3 in symmet- 
ric (22.5° — 22.5°) dd junctions as opposed to the usual 
estimate for a GBJ, D ~ 10~ 5 — 10~ 2 . Since usually 
I]} III oc D, a high-transmissivity GBJ is required in 
order to observe effects of the second harmonic. An esti- 
mate of the average transmissivity of our junctions would 
be p a bl/B,NA ~ 10~ 2 assuming I, the mean free path, 
to be equal to 10 nm and a resistivity in the a-b plane p a b 
equal to 10 -4 ilem. This is still too low to explain the 
strong 2 nd harmonic we observe. However, it is known 
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from, e.g., TEM-studies |19|. that the grain-boundary is 
far from uniform; the properties can significantly vary 
depending on the local properties of the interface, effects 
such as oxygen diffusion out of the GB etc., which are 
difficult to control. It is therefore reasonable to assume 
that there are many parallel transport channels through 
the GBj^,!^]. Channels with high transmissivity dom- 
inate the transport and might have D ~ 0.1 even though 
the average transmissivity is much lower. This is also 
consistent with the fact that most of our SQUIDs seem 
to be highly asymmetrical which is to be expected if the 
distribution of channels is random. The ratios of l\ and 
l]} can vary as much as ten times between two junctions 
in the same SQUID, even though the fluctuations of the 
total I c from sample to sample are much smaller. It is 
also clear from general considerations that a high value 
of I^ 1 excludes a high value of ij., since the 2 nd harmonic 
usually dominates if the odd harmonics of the supercur- 
rent are cancelled by symmetry |2flj . 

Recent studies of 0° — 45° GBJs have demonstrated 
that the SQUID-dynamics can be altered by the d-wave 
order parameter in YBCO |30| . It is however important 
to point out that our results do not directly relate to e.g. 
tetracrystal 7T-SQUID experiments; the latter crucially 
depend on having one 7r-junction with negative critical 
current, but still only the first harmonic present in I c (4>). 
Our SQUIDs have a conventional geometry, but uncon- 
ventional current-phase relations. 

One explanation for the pronounced effects of the 2 nd 
harmonic could be that relatively large sections of the 
interface are highly transparent and have a low degree of 
disorder. This in turn could be related to our fabrication 
scheme which seems to preserve the integrity of the bar- 
rier. This makes feasible their applicability in the quan- 
tum regime and supports our expectations that quantum 
coherence can be observed in this kind of structures. 

To summarize, we have observed very pronounced 2 nd 
harmonic in the current-phase relation of a 'conventional' 
YBCO dc-SQUID with 0° - 45° grain boundary junc- 
tions. It has strongly influenced the SQUID dynam- 
ics. All details of the SQUID behavior were explained 
within a simple model of a dd-junction with relatively 
high transparency. We believe that these effects are im- 
portant for better understanding of HTS Josephson junc- 
tion and SQUIDs. 
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